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Abstract 

We present direct measurement of Kirkwood-Rihaczek (KR) distribution for spatial properties 
of coherent hght beam in terms of position and momentum (angle) coordinates. We employ a 
two-local oscillator (LO) balanced heterodyne detection (BHD) to simultaneously extract distri- 
bution of transverse position and momentum of a light beam. The two-LO BHD could measure 
KR distribution for any complex wave field (including quantum mechanical wave function) without 
applying tomography methods (inverse Radon transformation). Transformation of KR distribu- 
tion to Wigner, Glauber Sudarshan P- and Husimi or Q- distributions in spatial coordinates are 
illustrated through experimental data. The direct measurement of KR distribution could provide 
local information of wave field, which is suitable for studying particle properties of a quantum 
system. While Wigner function is suitable for studying wave properties such as interference, and 
hence provides nonlocal information of the wave field. The method developed here can be used for 
exploring spatial quantum state for quantum mapping and computing, optical phase space imaging 
for biomedical applications. 

PACS numbers: 42.50.Dv, 42.65.Lm, 03.67.Hk 
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I. INTRODUCTION 



Optical phase-space tomography is an optical method for charac- 

terizing spatial properties of light fields/photons. The Kirkwood (K(x,p)) [9] and Wigner 
(W(x,p)) [l^ distribution functions were originally proposed in the studies of quantum 
statistics and thermodynamic for almost classical ensembles. The Kirkwood distribution 
('0(x)'0*(]9)e~^^^)has been rediscovered by Rihaczek for use in the theory of time- 
frequency analysis of classical signals. The Wigner function is more popular than the KR 
function because it has some unique properties such as negative values, real and symmetry. 
It can be obtained through tomograph methods and applicable to problems in phase space 

12| . The Wigner function was first intro- 



transport equation such as Liouville equation 
duced in optical light field by Bastiaans 
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lj| to analysis spatial properties of an optical 



Gaussian beam. In quantum optics, the quantum mechanical wave function cannot be mea- 
sured in principle. The wave function or quantum state of a physical system can be best 
represented by the Wigner function. Vogel and Risken [ist had theoretically proposed how 



quadrature amplitudes of nonclassica. 
by tomograph methods. Raymer 



li ght fields can be represented in Wigner distribution 



11, 
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has pioneered the tomography measure- 
ment of Wigner function for quadrature amplitudes of squeezed light, spatial properties and 
time-frequency properties of coherent light (coherent state with large mean photon number). 
In most quantum information experiment, a single spatial mode of electromagnetic field 



is generally used and considered [21 



221] . Quantum fiuctuations of light at diflFerent spatial 



points in the transverse plane of the light beam have to be taken into account. Kolobov 



and Sokolov 
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24. 



27l | have studied in detail the possibility of local squeezing in a 



spatial mode of light beam. Multimode spatial modes would cause information processing 
and computing errors. The property of multimode squeezed light which allows us to increase 
the sensitivity beyond the shot-noise limit, could create many interesting new applications in 
optical imaging, high-precision optical measurements, optical communications and optical 
information processing. However, one has to employ a dense array of photodetectors to 
observe multimode squeezed states. 

The phase space distributions from the Glauber Cahill s-parameterized class of 
quasi-distributions that contain the Wigner function, the Glauber Sudarshan P- 
representation {29^, and the Husimi or the Q-representation, have been widely used 
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as powerful phase space tools. However, most of tomograph methods are not involving di- 
rect measurement of these quasi-quantum distribution functions. These distributions are 
usually reconstructed from the raw data, where numerical and mathematical transforma- 
tion discrepancies may have made these distributions doubtfully representing the true states. 
The direct measurement of these distributions are desirable to represent the true quantum 
state. It is well understood that tomograph methods developed for these distributions can 
be used for any complex wave fields regardless of classical or quantum origin 3lj] . In this 
perspective, we use wave field to represent quantum mechanical wave function and coherent 
hght field. 

The KR distribution for any quantum state in a generalized form has been introduced [33|]. 
The hydrogen atom has been investigated using the Kirkwood— Rihaczek phase space rep- 



resentation 



32|. It has been pointed out that the KR distribution may not suitable for a 



"rotated state". The Wigner function for a rotated state is just simply the rotated Wigner 
function. This is the basis property in tomography method applied in reconstruction of 
Wigner function. However, this simple relation does not hold for KR function. Tomograph 
methods (Radon Transform) are not suitable for KR function. It is worth noting that posi- 
tion and momentum distributions of wave field can be obtained through the KR distribution 
such as J K{x^p)dp = and J K{x^p)dx = respectively. These properties are 

similar to Wigner function. Moreover, the complex conjugated KR can be directly measured 
through a two-LO technique [1], not like Wigner function which is reconstructed from ex- 
perimental data using quantum tomographic methods. Phase-space distribution functions 
especially Wigner function can unravel unique quantum properties such as entanglement of 
correlated systems 0, 35], negative parts of phase-space plots jsol, and the phase space 



sub-Planck structures of quantum interference |37l. I38l|. The KR is relatively unexplored for 
quantum fields. 

Spatial coherence of light sources is necessary for achieving efficient coupling into fiber 
systems for quantum communication, and also for biomedical imaging in image-guide inter- 
vention (sol. From a quantum-mechanical perspective, spatial degree of freedom of a photon 
is another optical quantum realization for encoding information such as spatial qubits (2l| . 
Full characterization of arbitrary, continuous spatial states of photons is important for un- 
derstanding the concept of the photon wave function [31]. Wigner function for an ensemble 
of identically prepared photons in the transverse spatial modes can be completely used to 
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characterize the transverse spatial state of the ensemble. The KR distributions are relatively 
unexplored for spatial properties of wave fields. In the second quantization of quantum me- 
chanics, electric field is written as operator in term of harmonic oscillator basis, that is, 
position and momentum operators. The spatial property of the electric field perpendicular 
to the propagation direction is described by mode functions. The mode function in x or p 
coordinate is a description of probability amplitude to find a photon at transverse position x 
or transverse momentum p, respectively. The field operator a{x) corresponding to the mode 
function will provide mean field and quantum noise Aa(x) in an optical detection 

scheme. In this work, distribution of x and p of a coherent wave field is measured by using 
a two-LO balanced heterodyne detection. Mean value measurement of the wave field 
at a configuration space will provide classical-like feature regardless the origin (classical 
or quantum light) of the measured field. Variance measurement of the wave field Aa(x) 
will provide quantum feature of the light field. For coherent state with low mean photon 
number or larere mean photon number, the variance for the quantum noise of wave field is 

n 

constant [40], that is, position or momentum independence. Therefore, the measurement 
method we have developed for coherent light is also applicable to coherent state, making 
them a useful testing ground for quantum imaging and mapping for information processing 
and quantum communication. Recently, EPR entanglement in spatial coordinate has been 
demonstrated by mixing an optical coherent light beam with squeezed light {4!]. However, 
their measurements had involved displacement and 'tilt' (momentum or angular) of a whole 
optical beam, not the transverse amplitude and phase structure of the optical beam. There 
is no approach in continuous variable quantum mechanics to measure the Wigner function 
or density matrix of spatial properties of EPR entangled beams. The optical technique 
and procedure presented in this paper could provide detail studies of phase-space physics in 
quantum metrology through sub-Planck phase-space structures in the Wigner function, and 
also discrete Wigner function for quantum mapping. 

Phase-space distributions are used to represent quantum-mechanical operators in explor- 
ing phase-space quantum effects and quantum-classical correspondence. Quantum algo- 
rithms for measuring KR and Wigner distributions have been developed Q • Discrete 
phase space distribution has been suggested to show potential advantages in quantum com- 
puting especially quantum mapping. Quantum algorithm can be simply thought of as a 
quantum map acting in a Hilbert space of finite dimensionality. Specifically, algorithms be- 
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come interesting in the large N limit i.e., when operating on many qubits!. For a quantum 
map, this is the semiclassical limit where regularities may arise in connection with its clas- 
sical behavior. These semiclassical properties may provide hints to develop new algorithms 
and ideas for novel physics simulations. 

In this paper, we demonstrate direct measurement of the KR distribution for a wave 
field with Gaussian mode function of TEMoo- Linear transformation from KR distribution 
to Wigner, P and Q-distribution are also plotted to show the fundamental differences in 
theirs respective. Second, a superposition of two spatially separated coherent light beams 
is used for discussing the phase-space interferences in these distributions. We show that 
complex conjugated KR distributions for spatial properties of wave fields can be determined 
by use of a novel two-LO balanced heterodyne detection scheme f])]. The technique measures 
/C*(x,]9) =< £''{x)£{p) > e^^^, which can be written as SR + iSj. A lock-in-amplifier is used 
to measure the Sr and Sj with respect to relative phase setting of a reference signal at 0° and 
90°, respectively. By changing the lock-in-amplifier reference phases such as 0° and —90°, the 
system will measure ]C{x^p). However, we keep the lock-in-amplifier with reference phases at 
0° and 90° for all measurements in this paper. There is no different in physics implied by wave 
field in KR or complex conjugated KR representation. The two-LO heterodyne technique 
was originally designed for biomedical imaging, that is for optical phase space coherence 
tomography of the light transmitted through or refiected from biological tissue. Now, we 
use this measurement to explore the properties of one particle wave mechanics or wave 
field through KR, Wigner, P- and Q- distributions. The technique can be used to measure 
any complex spatial wave fronts such as divergence and convergence, and phase-conjugated 
properties of wave field . The /C*(x,]9) can be easily transformed to Wigner function 
by using a linear transformation where the Radon transform is not required. The advantage 
of KR is it contains local information of the wave field. If there is no wave field presents at 
a configuration space (x,p), then there will be no distribution at the (x,p). It serves better 
for optical imaging in biomedical application such as to characterize the cell structure. We 
will illustrate this spatial property of KR distribution and compare it with Wigner function, 
P- and Q- distributions. In the two-LO balanced heterodyne detection, we use a local 
oscillator (LO) field comprising a coherent superposition of a tightly focused LO Gaussian 
beam of TEMoo and a highly collimated LO Gaussian beam of TEMoo- This scheme permits 
independent control of the x and p resolution, permitting concurrent localization of x and 
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p with a variance product that surpasses the minimum uncertainty hmit associated with 
Fourier-transform pairs. Quantum mechanics did not allow simultaneously measurement of 
X and p of a wave field. However, simultaneous measurement in distribution of x and p for 
a wave field is allowed. 

II. CHARACTERISTIC FUNCTION APPROACH 

In this section, the characteristic function method will be used to transform the charac- 
teristic function of KR distribution to Wigner, P- and Q- distributions. The characteristic 
functions in term of harmonic oscillator basis will be used through out this paper because 
it is more revelent to spatial properties of coherent wave fields. Since we directly measure 
the complex conjugated of KR distribution, we will discuss how the measured results can be 
used to obtain Wigner, P- and Q- distributions. We start with the characteristic function 
for the complex conjugated of the KR distribution in term of coherent state representation, 
as given by, 

Mkr{P, n = Tr{pe-\^'-^'"\l'e'^^'^'P^^). (1) 

where (/?, /?*) are the Fourier conjugate pairs for (a, a*)the eigenvalues of a and a*, respec- 
tively. The interesting feature of this characteristic function in the Fourier plane Mkr{P, /?*) 
is the trace of displacement operator T>{f3) followed by squeezing operator <S(l/2). The ab- 
solute of 1/?^ — /?*^ 1/4 just to avoid the confusion of other forms of definition that are, f3 if3 
and /?* — z/?*. The complex conjugated of KR distribution in the (a, a*) plane can be 
obtained through, 

)C*{a,a*) = d'f3e-'^"*-'^*"MKRiP,n- (2) 

The /C*(a,a*) can be transformed to Wigner, P- and Q- distributions in term of (a, a*) 
through the relationships of the characteristic functions, such that, 

Mkr{P, n = e-\P'-P*'\/'Mw{(3, n (3) 

where, 

MwiP, n = e-\^\"'''Mp{(3, P*) = el'^l'/^A^Q. (4) 
The characteristic function for Wigner function is given by, 

MwiP, (3*) = rr(pe'^"^+*^*"). (5) 
6 



The P- and Q- representations are related to the characteristic function of Wigner function 
through normal ordering and anti-normal ordering of (a, a*). 

In our experiment, we measure spatial properties of a wave field in term of position and 
momentum coordinates (x,p). The ]Q{x^p) can be obtained from Eq. ([T]) and Eq. ([2j) using 
the variables, 

^=-^(a + zr);/?* = -^(a-zr) 

a = -^(x + zp); a* = -^(x - zp). (6) 

Then, the following terms, 

^—if3(j*—i(3*cr ^—iax—irp ^ij^ 

are obtained. The ]Q{x^p) can be written as, 

^.(.,rt^/d.*e— ^T.(pe-'-e-). (8) 

By using the identity J dx\x){x\ = 1 and p = we obtain, 

/C*(x,p) = / e-'^Pij{x + T)dr. (9) 



By changing the variable ^ = x + r, rf^ = rfr, we obtain, 

/C*(x,j9) = -V^*(x)V^(p)e^^^. (10) 

TT 

which is the complex conjugated of KR distribution. In order to write the characteristic 
functions of Wigner, P- and Q-distributions in term of spatial properties of coherent light 
beam such as beam waist a, position and momentum coordinates (x,p), we use the variables, 

so that the characteristic functions for KR, Wigner, P- and Q- distributions are related to 
each other as given by, 

Mkr{x,p) = e-'-'p'/''Mw{x',p') (11) 



and, 

Mw{x,p) = e-^("'^^'+"^'/"')>^p(x^y) = ei("'^^'+^^'/"')>^g(x^pO• (12) 

Since our experiment measures the /C*(x,]9), its characteristic function is obtained through 
transformation as given, 

Mkr{x'p') = j dxdp]C\x,p)e'''P'^'^''' . (13) 
Then, the Wigner function is obtained through, 

or in the simphfying form as in the Ref i.e. 

W{x,p) = ^ J rfx>'e-'^(" -")(^-^)/C*(x,p). (15) 
Similarly, the P- and Q- distributions can be obtained through, 

P{x,p) = J e'^^'^'/^MKR^x ,p')e^^^"^'"^'''"'^"^^ (16) 

and, 

respectively. 

The similarity of KR and Wigner function is when they are integrated over mo- 
mentum/position, the two functions will provide the same result for the probability in 
position/momentum i.e. j ]Q{x^p)dx = JW{x^p)dx = \ip{p)\'^ and j ]Q{x^p)dp = 
J W{x^p)dp = |'0(x)p. It should be noted that because of the uncertainty principle, nei- 
ther function has a physical meaning until it is integrated over either momentum space or 
configuration space. 

III. EXPERIMENT DETAILS 

Since tomography methods developed for spatial properties of photon wave function can 
be applied to classical field (coherent state with large mean photon number), we use electric 
field notation S{x) to represent the wave field ip{x) in the following Sections. 
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signal beam 



FIG. 1: A balanced heterodyne optical imaging scheme. 

We use a balanced heterodyne detection scheme as shown in Fig. [H The beat amplitude 
Vb is determined by the spatial overlap of the local oscillator (LO) and signal (S) fields in 
the plane of the detector at Z = Zjj as, 

Vb = J dx'£lo{x\zD)Ss{x',ZD) (18) 

where x' denotes the transverse position in the detector plane. When the LO beam is moved 
by a distance the Vb becomes, 



Vb(4) = j dx'8lQ{x' - 4, zd)Ss{x', zd)- 



(19) 



The fields in the detector plane are related to the fields in the source planes [z = 0) of lenses 
LI and L2, which have equal focal lengths f = 6 cm. The LO and signal fields at 2; = after 
the lenses, LI and L2, are given by 

k 

^Loi^i -do:,z = 0) = exp[-i—xl]£Lo{xi -d^,z = 0) 



k 

£'s{x2, z = 0) = exp[-i—xl]£s{x2, z = 0) 



When the lens L2 is scanned by a distance dp, the signal field (EOl) is altered as 

k 



(20) 



£'s{x2, z = 0) = exp[-i^(x2 - dpf]£s{x2, z = 0) 
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The fields propagating a distance d = f to the planes of the detectors can be obtained by 
using Fresnel's diffraction integrals as, 



/ k f k 
^Lo{x' - d^, zd) = ]J-^ J dxi exp[i—{xi - x'Y] 

k 

X exp[-i—xl] Slo{xi -d^,z = 0) 
Ss{x',zd) = )J~^J J dx2 exp[i-^{x2 - x'f] 



k 

X exp[-i—{x2-dpf]£s{x2,z = 0) 



(21) 



By substituting the above equations into Eq. (fT9l) . the quadratic phases in cancel and 
the quadratic phases that depend on 2 cancel in these expressions because the detector 
plane is in the focal plane of the lenses, LI and L2. One obtains 

k k f k 

Vsidx^dp) = y exp[i— rfj] j dx2 exp[-ijX2dp] Ss{x2, z = 0) 

X J dxi £lo{xi -d^,z = 0)5{xi - X2) (22) 

Integrating over Xi and by replacing X2 by x and dropping the term z = 0, the mean square 
beat amplitude is then given by 

2 



\VB{d^,dp)\^ oc 



J dxSlQ{x - dx)Ss{x) exp(-i/cyx) 



or 



/d 
dx£lo{x - d^)£s{x) exp(-z/cyx) 

X j dx' £lo{x' — dx)£s{^') exp(i/cyx^). 



(23) 



This can be rewritten using the variable transformations x = Xo + r}/2] x' = Xq — r]/2 
where the Jacobian of this transformation is 1. Then, by using the definition of the Wigner 
distribution, 

W{x,p) = j ^exp{i€p){£*{x + €/2)£{x-€/2)) (24) 
and its inverse transform is given by, 



£*{x + e/2)£{x - e/2) = j dp exp{-iep)W{x,p), 



(25) 
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the beat signal in Eq. (031) becomes 

|Vb(4,c?p)P oc j j dr]Elo{xo + r]/2- d^)£Lo{xo- ri/2- d^) 

/d 
dp exp{—ik^rj) exp{—irip)Ws{x,p). (26) 

Since the Wigner distribution of the LO field is 

Wlo{xo - d^,p + k^) = J exp[-iri{p + k^)] 

X SloiXo + V/'^-d:c)£Lo{Xo-V/'^-d:c) (27) 

then the LO fields in Eq. f l26l ) can be replaced by the Wigner function in Eq. f l27l ). Finally, 
the mean square heterodyne beat signal \Vb\'^ can now be written as 

/k 
dxdpWLo{x - d^,p+ —)Ws{x,p) (28) 

where Ws{x^p) is the Wigner distribution of the signal field in the plane of L2 (z =0) and 
Wlo{XiP) is the LO Wigner distribution in the plane of LI (z = 0). We include a detail 
description of two- window heterodyne measurement of KR distribution as shown in Fig. [2l 
The variables dx and dp respectively indicate the positions of a mirror Ml and a lens L2 
as in Fig. [2l Eq. f l28l ) shows that the mean-square beat signal yields a phase-space contour 
plot of Ws{x^p) with phase space resolution determined by Wlo{XiP)- By using a two-LO 
heterodyne detection scheme as discussed below, the \ VB{dx^ is found to be proportional 
to /C*(x,p). 

To obtain independent control of the x and p resolution in heterodyne measurement, we 
employ a slowly varying LO field containing a focused and a collimated field with a well 
defined relative phase 9 

£lo{oc) = fo[exp(- — ) + a exp(-^) exp{ie)] (29) 

Here a is chosen to be small compared with the distance scales of interest and 1/A is 
chosen to be small compared with the momentum scales of interest in the signal field. The 
schematic picture of the overlapping of the LO beam with the spatial width a and another 
LO beam with the spatial width A is shown in Fig. [3l One can see that the overlapping 
area is determined by the position and momentum resolutions for the LO fields in Eq. f l29l ) . 
The focussed LO gaussian beam extracts the position information of the signal field and 
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FIG. 2: Experimental setup for measuring the KR distribution using two-LO balanced heterodyne 
detection. X^; squarer, D; photodiode detector, BS; beam splitter. 

the coUimated LO gaussian beam extracts momentum information of the signal field. The 
Wigner function for the LO field is obtained by substituting Eq. f l29l ) into Eq. f l24l ). We take 
J^^a^ . In this case the phase- (0-) dependent part of the Wigner function for the LO takes 
the form 



- 2aV) cos(2xp + Q) 



cos{2xp 



(30) 



where the last form assumes that the range of the momentum and position integration in 
relation (l28l) is limited by the signal field. 

The measurement of phase space distributions is accomplished by translation of optical 
elements. These elements are all mounted on translation stages driven by computer con- 
trolled linear actuators. The system scans the LO position over a distance dx = ±1 cm by 
translating mirror Ml in the LO path. The LO momentum is scanned over ±0.3 fc, where 
A; = 27r/A is an optical wave vector, by translation of the signal-beam input lens L2 (focal 
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FIG. 3: The overlapping area is the position and momentum resolutions of the combined LO 
beams. The uncertainty is proportional to -| less than 1. 

length / = 6 cm) by a distance dp. 

In the experiments, as illustrated in Fig. [21 the LO beam is obtained by combination 
of two fields that differ in frequency by 5 kHz, so that 9 =(27r x 5kHz)t. Lens L3 focuses 
beam LOl to a waist of width a, and lenses L4 and L5 expand beam L02 to width A. We 
combine these two components at beam splitter BS3 to obtain an LO field of the form given 
in Eq. f l29l ). We monitor one output of the beam splitter with detector to phase lock 
the 5 kHz beat signal to the reference channel of the lock-in amplifier. Each component of 
the LO beam is shaped so that it is at a beam waist at the input plane of the heterodyne 
imaging system (lens LI). The focussed gaussian LO beam is frequency shifted at 110 MHz 
and the coUimated gaussian LO beam is frequency shifted at 110 MHz plus 5 kHz. These 
two LO beams are overlapped with each other and phase locked at 5 kHz. The signal beam 
is frequency shifted at 120 MHz. Two imaging lenses LI and L2 are used to overlap the dual 
LO beam with the signal beam at two detectors. The interference beat signal between the 
signal beam and the dual LO beam is obtained at detectors 1 and 2 and consists of 10 MHz 
and 10 MHz plus 5 kHz components. These signals are sent to a spectrum analyzer. The 
spectrum analyzer bandwidth, 100 kHz, is chosen to be large compared with 5 kHz difference 
frequency. The output of the analyzer is then squared by using an analog multiplier. The 
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mean square signal has components at 5 kHz. A lock-in amplifier is used to measure the in- 
and out-of phase components of the multiplier output at 5 kHz. The lock-in outputs for in- 
and out-of-phase quadratures then directly determine the real and the imaginary parts of 
the /C*(x,]9) function 

|^B(4,rfp)r = /C*(^o,Po) = J ^^^^ exp[2i{x' -Xo) X {p' -p,)]Ws{x',p') 

= {S^{xo)S{po)) exp{iXoPo) = Sr + iSi, (31) 

where the Wlo{x^p) in Eq. (I30l) is replaced by e^(^^^+^) in Eq. (I28l) yielding the in-phase and 
out-of-phase contributions from 9. Here Xo = dx is the center position of the LO fields and po 
= —kdp/f is the center momentum. The Sr and Sj components are related to heterodyne 
beat signal of E^q^Es and EIq2Es or intensity correlation of E^q^Es and EIq2Es. In the 
balanced heterodyne detection, the output voltage VB{t)=Vi{t) — V2{t) before being fed into 
the spectrum analyzer is given by, 

Veit) = 2ElELoie~'^^'^'''~^'^' + 2ElEL02e~'^^'^''^~^'^^ + cc (32) 
In spectrum analyzer, the power spectrum of the Vb is measured as 

^e'^^{VB{t)VB{t + T)). (33) 

As mentioned previously, it is squared by using a squarer to recover the beat signal |V0(f^)p. 
From Eq. f l33l ). the power spectrum can be calculated by keeping the slowly varying term 
(^LOi — ^L02) in time t and other terms that depend on r, that is. 



{VB{t)VB{t + T)) = ^ r'\tVB{t)VB{t + T 
^ J-T/2 



-r/2 

+ (: the negative frequency contributions from the above terms:) 

(34) 

Here, Qs=120 MHz, ^^01=110 MHz + 5 kHz and ^^02=110 MHz. Now, by substituting 
Eq. ( IMj) into Eq. (jSHl) to obtain the power spectrum for the beat Vb and setting the analyzer 

14 



at 10 MHz with the bandwidth of 100 kHz, the |K>(t)P at 5 kHz after the recovery by the 
squarer is 



\Vo{t)\' (X E'sELOiEsElo,e 



i{^L01-^L02)t 



+ CC 



(35) 



Here Qloi — ^lo2 = 5 kHz. The in-phase and out-of-phase components of the 
at 5 kHz correspond to Sr and Sj in Eq. (IHTI) . Note that E^Elqi is integrated over the 
transverse plane as is EsEIq2. It is worth noting that the component E^E^oi of Eq. (I35l) 
is corresponding to the measurement of the position distribution in K*(xo,Po) of Eq. (ISB 
by the tightly focussed LOl beam. Similarly, the component EsEIq2 of Eq. f l35l ) indicates 
the measurement of the momentum distribution in X*(xo,Po) of Eq. (IHTI) by the collimated 
L02 beam. 

As the position of mirror Ml is scanned a distance dx^ the optical path lengths of the LO 
fields change. For the current experiments, the HeNe laser is a source, the change in path 
lengths is small compared with the Rayleigh length and the coherence length of the beams, 
so translating Ml simply changes the center position of the LO fields. 

IV. RESULTS 

A. Measurement of an Optical Gaussian beam 

1. Gaussian Beam 

As an initial demonstration of the capability of this system, we measure the /C*(x,]9) 
function for an ordinary gaussian beam. A one dimensional wave field for an Gaussian beam 
of TEMqo and radii of curvature, R, is given by. 



where the x is the transverse position and a is the waist of signal beam. The Fourier 
transform of S{x) is. 



£{x) oc exp( 




(X A + iB 



(36) 



£{p) oc exp( 



oc C + iV 



(37) 
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The complex conjugated KR distribution for this wave field can be written as, 

/C*(x,j9) = {AC + BV) co8{xp) + {BC - AV) 8m{xp) 

+ i{{AV - BC) co8{xp) + {AC - BV)sin{xp)) (38) 

The characteristic function, Aijcn{XiP)i is obtained from Eq. (fTTI) . where the characteristic 
function for Wigner function for this wave field is given by. 



I 



Then, the Wigner, P- and Q-distributions are obtained from Eq. ([HI), Eq. (fT6l) and Eq. (fTTI) . 
respectively. The Wigner function for the wave field is given by, 

W{x,p) = iexp(-4 - tH^+P?) (40) 
The P- and Q-distributions for the wave field are given in the integral form as, 

P{x,p) = / exp( — )exp{—ixp — ipx )ax dp 

J AR 2R 

Q{x,p) = / exp(-— — ^^)exp{-txp -ipx)dxdp 

(41) 

For simplicity, the signal Gaussian beam is shaped by a telescope so that its waist coincides 
with input plane L2 of the heterodyne imaging system. For a gaussian beam at its waist, 
i? = oc, Eq. (I38l) gives the complex conjugated KR distribution as, 

/C* (x, p) = AC cos(xp) + iAC sin(xp) (42) 

2 2 2 

where A = exp(— ^), C = exp(— ^^), and (7=0.85 mm is the 1/e-intensity width. The 
]Q{x^p) function for the signal field is measured by use of the dual LO beam of the form 
given by Eq. f l29l ) with a =81 /xm, A = 2.6 mm, and a =1. The measurement result for a 
gaussian beam is shown in Fig. [H The top row is our experimental results and the bottom 
row is a theoretical prediction obtained by using Eq. f l42l ). The real and the imaginary parts 
of the detected signal, Eq. f lHTl ). are shown in Fig. [4](a) and (b). Our observation is similar 
to the theoretical prediction by Wodkiewicz js^ for a coherent state. 
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FIG. 4: The Wigner function for a gaussian beam. Top row is experimental results and bottom 
row is theoretical prediction for a gaussian beam, (a) is in-phase component of /C*(xo,Po) and (b) 
is out-of-phase component of /C*(xo,Po) 



Arbitrary (a) 

0.005 I , , , 



0.004 - 
0.003 - 
0.002 - 
0.001 - 



-3-2-10 1 2 3 

Position (mm) 



Arbiti 

0.004 



-0.002 
-0.004 



(b) 



Position (mm) 



FIG. 5: (a) Real part and (b) imaginary part of position distribution of a Gaussian field. Imaginary 
part is around zero. 

Position (momentum) distributions of this field can be obtained through the summation 
of momentum (position) coordinate of real and the imaginary parts of the measured ]Q{x^p). 
The position and momentum distributions are shown in Fig. [5] and Fig.[6l respectively. The 
imaginary part of position and momentum distribution are around zero as theoretically 
predicted by J IC^{x^p)dp = \£{x)\'^ and J JC^{x^p)dx = \£{p)\'^^ respectively, which are 
the real physical quantities (no complex values). From here, the position and momentum 
distributions are fitted with Gaussian function. We obtain the beam waist of a = 0.86 
from position distribution and cr=0.87 from momentum distribution. Both results are in 
excellent agreement with the measured width cr=0.85 mm obtained by use of a diode array, 
demonstrating that high position and momentum resolution can be jointly obtained. 

Wigner distribution is obtained by using simple linear transformation as discussion in 
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FIG. 6: a) Real part and (b) imaginary part of momentum distribution of a Gaussian field. Imag- 
inary part is around zero. 

Eq. ( fT4l ) or Eq. f fTSl ) as shown in Fig. [71(a). We fit the width of the measured in-phase 
signal in position for p = to obtain a spatial width a =0.87 mm, whereas the corresponding 
momentum distribution for x = yields a =0.83 mm. 

The characteristic function A^xi?(^^pO obtained by using numerical integration of 
the measured IC^{x.p) as in Eq. (fT3l) . The P- and Q- distributions for this signal gaussian 
beam are then obtained through the Eq. f fT6l ) and Eq. f fTTI ). as shown in Fig. [71(b) and 
Fig. [71 (c), respectively. The P-distribution has a narrower peak in phase-space compared 
to other distributions. This is predicted for the signal beam with i? = cxd in P-distribution 
of Eq. dm), which we should have 5{x)5{p) in phase-space. The Q-distribution has a broad 
peak in phase-space compared to other distributions. The Q-distribution for this signal 
beam can be evaluated from Eq. (HB with i? = oc, as given by, 

2^2 ^2 

Q{x,p) ocexp( — ). (43) 

The position width of the Q-distribution at p=0 is about ^/2 larger than the position width 
of Wigner distribution at p=0 as from Eq. f l40l ). Hence, the beam waist for the signal beam 
obtained from the Q-distribution is ^/2 larger than the exact value. 

2. Measurement of superposition of two slightly displaced coherent beams 

As a fundamental feature in the process of quantum measurement, we cannot observe 
physical properties of a quantum objects directly because the overall backaction of any 
observation cannot be made less than Planck's constant h. Instead, we observe the wave or 
the particle aspects of the physical objects. The ]Q{x^p) distribution for a coherent field is 
more likely representing particle picture of the field because it contains local information of 
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FIG. 7: (a) The reconstructed Wigner distribution, (b) Husimi or Q-distribution, and (c) Glauber 
Sudarshan P-distribution for an optical Gaussian beam. 2D-plot (left) ;3D-plot (right). 

the coherent field. While the Wigner function is more likely representing wave behavior of 
the field because it exhibits phase-space interference. Based on these properties, the ]Q{x^p) 
and Wigner distributions are very useful to characterize a wave field or an physical object 
through phase-space imaging in many applications such as quantum imaging, metrology, 
and biomedical imaging. To illustrate the particle picture of ]0{x^p) and the wave picture 
of Wigner function, we use the same signal gaussian beam obscured by a wire with diameter 
of 1 mm. Then, the electric field ^s(x) as a function of position is shown in Fig. [8l It 
didn't involve convolution integration because the wire is placed close to the imaging lens 
L2. It shows the superposition of two slightly displaced coherent beams. It is analogous 
to a Schrodinger cat state. In this case, the slowly varying field is gaussian as before but 
multiplied by a slit function that sets the field equal to zero for \x\ < 0.5 mm. Fig. [9] shows 
contour plots of the real and the imaginary parts of the detected signal of /C*(x,]9). 

The top row is our experimental results. The bottom row is theoretical prediction. The 
theoretical plots are first obtained by numerically generating the signal function £s{x) and its 
Fourier transform Ss{p)- The real and imaginary parts of /C*(x,]9) in Eq. f fTOl ) are then theo- 
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FIG. 8: The wave field £{x) of a Gaussian beam obscured by a wire. 
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FIG. 9: Top row shows experimental results and bottom row shows theoretical predictions for a 
gaussian beam blocked by a wire, (a) in-phase component of the /C*(xo,Po) and (b) out-of-phase 
component of the /C*(xo,Po). 

retically plotted. The real and imaginary parts are not showing phase-space interferences of 
two slightly displaced coherent beams. These plots show local information of the signal field, 
i.e., zero(nonzero) field is corresponding to zero(nonzero) 10{x^p) phase-space distribution. 
This locality property exhibits particle picture if an atomic wave function/single photon 
function is used. The position and momentum distributions are then retrieved J ]0{x^p)dp 
and J 10{x^p)dx as shown in Fig. [TOland Fig.dU respectively. The momentum distribution 
contains the interference features of two spatially separated wave packets of £s{x) as shown 
in Fig. [m The imaginary part of these distributions are around zero as expected. 

The Wigner function as shown in Fig. [12] is reconstructed by using the linear transforma- 
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FIG. 10: (a) The real part, (b) the imaginary part of the position distribution of a wire function 
obtained by integrating the measured KR distribution over momentum. Imaginary part is around 
zero. 
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FIG. 11: (a) The real part, (b) the imaginary part of the momentum distribution of a wire function 
obtained by integrating the measured KR distribution over position. Imaginary part is around zero. 

tion of the measured /C*(x,]9). The coherence between these two wave packets in the signal 
field leads to an phase-space interference pattern in the momentum distribution. The sig- 
nature of this coherence in the Wigner distribution is the oscillating positive and negative 
values between the main lobes. An interesting feature of this Wigner distribution is the 
oscillation in momentum at the position x =0 of the wire. We observe the negative values 
which is analogous to quantum interference in phase space. This feature can be seen in 
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FIG. 12: The reconstructed Wigner distribution for a wire function (a) two-, (c) three- dimensional 
plot. 
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Fig. [121 which the reconstructed Wigner function is shown as a three dimensional plot. 
The negative values highlight the impossibility of a particle simultaneously having a precise 
position and momentum. It also makes sure that the sum over momentum along x = in 
the reconstructed Wigner distribution has zero intensity at the center. The negative and 
positive parts of the Wigner phase space distribution are important features to obtain full 
information about the field. Our observation of the negative values of Wigner function for 
this coherent field did not claim that the negative values exhibited by a quantum field is 
not a quantum feature. We believe that classical or quantum features of an experiment are 
based on classical or quantum field involved in the experiment. 

From the obtained Wigner function, one can obtain the position and momen- 
tum distribution of the Ss{x) by using the formulas \£s{p)\'^=JT^Ws{x^p)dx and 
\£s{x)\'^=J^^Ws{x^p) dp as shown in Fig. [T3l(b) and (d) respectively. The resolution of 
these plots is better than Fig. [10] and Fig. [TT] because we numerically generate more data 
points from the reconstructed Wigner function. The measurements are in agreement with 
the theoretical predictions as shown in Fig. [T3l(a) and (c) respectively. 

As discussed before, we first obtain the characteristic function M.kr{x\p') and then the 
P- and Q- distributions for this signal field. These distributions are plotted in Fig. [HI (a) 
and (b), respectively. The Q-distribution exhibits the broadening or low resolution of phase- 
space features compared to Wigner function. While the P-distribution is not a well behaved 
function for this field as it is expected for P-distribution. 



V. DISCUSSION AND CONCLUSIONS 



Spatial properties of photon wave mechanics approach [31] for a single photon have been 
studied in detail. The approach can be applied to both single photon state and coherent 
field. As a consequence, the coherent field is the best testing ground for developing tomog- 
raphy method for quantum information processing. In the similar eflForts, coherent fields 
have played an important role in quantum communication and computing such as search 



algorithm |44l, l45l. l46l . l47[ | and factorization of numbers [48|] . Optical wave mechanics imple- 
mentations [491, |50|] of entanglement and superposition with coherent fields (coherent state 



with large photon numbers) have been demonstrated. This implementation has been used 
to study entanglement swapping and tests of non-locality. However, photon wave mechanic 
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FIG. 13: The momentum and position distribution for the Schrodinger cat state, (a) and (c) are 
theoretical prediction of momentum and position distributions of the cat state, (b) and (d) are 
the corresponding experimental results of momentum and position distribution of the cat state 
obtained by integrating the measured Wigner distribution over p and x respectively. 

approach for two- and multi-photon spatial qubits haven't been relatively explored. The 
two-LO technique developed in this paper will provide another tomography tools to explore 
spatial qubits because it could provide the particle picture through KR distribution and 
wave picture through Wigner distribution. 

We would like to discuss KR and Wigner distribution separately because the particle 
and wave picture of wave field can provide independent useful information for quantum and 



39]. 



coherent information processing including biomedical imaging 

Particle picture or local information of a wave field represented by KR distribution has 
advantages in positioning or angle (momentum) resolving to extract local activities of an 
target. The method can directly locate the field or structure of an object without applying 
any raw data transformation. This will be particularly useful in cell tissue characterization 
such as prostate cancer cell detection. The KR distribution can be used to study Goos- 
Hanchen (GH) shifts [SJj] occur in near field optics and photonic waveguide. GH shifts in 
position and momentum can be used to identify the loss due to photonic crystal waveguide 
fabrication. 
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FIG. 14: The reconstructed (a) P- and (b) Q-distributions for a slightly displaced coherent beams. 
2D-plot(left); 3D-plot (right). 

Wave picture or nonlocal information of a wave field is best represented by Wigner func- 
tion. We observe the non-positive properties of the Wigner function for a superposition of 
two spatially separated gaussian field £{x) analogous to a Schrodinger cat state in spatial 
coordinate. We have demonstrated the similarities in phase space interference between spa- 
tial properties of quantum and coherent fields via the measurement of the Wigner function. 
We also show that an interesting analogy exists between our choice of LO field and that em- 
ployed in the quantum-teleportation experiment |52|. In the x-p representation, the small 
and the large beams of our two-LO can be viewed as superposition of the position (in-phase) 
and the momentum (out-of-phase) squeezed fields. A £{x) spatial gaussian field of TEMqq 
is the lowest mode and is similar to a coherent state in phase-space picture. Gaussian beams 
of smaller (larger) size than the lowest mode correspond to position (momentum) squeezed 
states. The two-LO technique here can only be used if we know the nominal size of the signal 
beam. The focussed and coUimated LO beams must be chosen to achieve sufficient x— and 
p— resolution for the given signal beam. These are the similar problems encountered in our 



In a complex multi-particle system or large N-biological system, we believe KR and 
Wigner distributions are very useful to study the local and nonlocal information of a macro- 
scopic wave field such as the mechanisms of decoherence due to neighbor particles, quantum 
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mapping of N-particle system or semiclassical system, and macroscopic entanglement be- 
tween two macroscopic mirrors. 

In conclusion, we have demonstrated the direct measurement of KR distribution using 
two-LO balanced heterodyne detection technique. The characteristic function of KR is 
related to Wigner, P and Q- distributions. Then, the Wigner, P and Q- distribution are 
plotted by using raw data from the KR distribution. The physical properties of a wave 
field such as local and nonlocal phase space information are illustrated through KR and 
Wigner functions, respectively. This two-LO technique can be used in information processing 
including quantum information for quantum mapping and optical imaging for biomedical 
applications. 
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